Classical Dynamics of Quantum Entanglement 



O 

(N 

Oh' 



Giulio Casati, 1,2 Italo Guarneri, 3 ' 4 and Jose Reslen 5 

1 CNISM, CNR-INFM & Center for Nonlinear and Complex Systems, 
Universita degli Studi dell'Insubria, Via Valleggio 11, 22100 Como, Italy 
2 Istituto Nazionale di Fisica Nucleare, Sezione di Milano, via Celoria 16, 1-20133 Milano, Italy 
Center for Nonlinear and Complex Systems, Universita degli Studi dell'Insubria, Via Valleggio 11, 22100 Como, Italy 
Jf Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, via Bassi 6, 1-21100 Pavia, Italy 
J Department of Physics, National University of Singapore, Singapore 117543 
(Dated: September 19, 2011) 

We numerically analyze the dynamical generation of quantum entanglement in a system of 2 
interacting particles, started in a coherent separable state, for decreasing values of K. As h — > the 
entanglement entropy, computed at any finite time, converges to a finite nonzero value. The limit 
law that rules the time dependence of entropy is well reproduced by purely classical computations. 
Its general features may then be explained by simple classical arguments, which expose the different 
ways entanglement is generated in systems which are classically chaotic or regular. 
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The development of quantum chaology has shown 
that quantum dynamics is deeply influenced by classi- 
cal chaotic motion. It was then natural to inquire, how 
is the dynamical generation of quantum entanglement 
affected by the chaotic or regular nature of the underly- 
ing classical dynamics [lrH, 0- 11 1 ; an especially intrigu- 



ing aspect of this question being that quantum entangle- 
ment does not have a classical counterpart. Since the 
early paper by Furuya et al. [l[, most investigations 
about this subject have supported the view that reg- 
ular classical dynamics leads to slower dynamical gen- 
eration of entanglement (DGE)Q. The opposite point 
of view, that DGE is scarcely affected by the nature 
of classical dynamics, has also been surmised Q; nev- 
ertheless there are convincing reasons to believe that 
the way quantum entanglement is dynamically gener- 
ated does indeed have a classical counterpart 0, @, 0|. 
Notably, DGE may be studied by resorting to a phase- 
space formulation of quantum mechanics [4j; and then 
the purity of reduced states, which is a measure of (dis- 
entanglement, can be directly obtained from the phase- 
space functions which represent such quantum states in 
the classical phase spaces. Letting such functions evolve 
according to the classical Liouville dynamics a classical 
counterpart to DGE is obtained In a similar line of 
thought, quasi-classical analytical methods have been im- 
plemented to investigate the decay of purity of reduced 
states in integrable systems 0, Q. 

In this Letter we study the h — > limit of DGE. We thus 
demonstrate that DGE has a direct classical analog, in 
the following sense. We start a system of two interacting 
particles in a separable coherent state, and numerically 
compute the growth in time t of their entanglement , as 
measured by the von Neumann entropy S(t, h) of their 
reduced states. In the limit H — > 0, the evolving quan- 
tum state tends to a classical pure state (i.e. a Dirac 
(^-function in phase space, supported in a point of the 



classical trajectory), and yet its entanglement entropy 
docs not tend to zeroQ; instead it converges to a finite 
value, and we find that this value can be approximated 
using Hamilton's equations. 

We consider a Toda model of two interacting particles 
moving in a line under the following Hamiltonian: 
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Relevant parameters are the energy E, and the mass ratio 
m\jm<x. For mi = m.2 the model is completely integrable 
no matter how large is the energy (Fig. la), while for any 
other mass ratio (except limit values 0,oo) the dynamics 
exhibits a transition from integrability to chaos as the 
energy E is increased [l2|. In particular for m%jm\ = 
0.54 motion is predominantly chaotic already at energy 
E = 7 (Fig.lb). 

To compute the quantum evolution we have first an- 
alytically calculated matrix elements of the quantum 
Hamiltonian (which is obtained from (fTJ) by canonical 
quantization) over the energy eigenbasis of two uncou- 
pled harmonic oscillators of unit mass, and then we have 
diagonalized the matrix thus obtained. With the cho- 
sen matrix size, the computed eigenbasis was verified to 
support the chosen initial states within negligible errors. 
Such initial states were separable coherent states, cen- 
tered at 91 = Q2 = and Pi = P2 = ^f~E in the regu- 
lar case (mi = 1712 = 1), and q\ = (72 = 0, p\ = \[E. 
P2 = —\pm2E in the chaotic case (mi = l,m2 = 0.54). 
The latter initial condition is inside the chaotic sea. As 
time goes on the particles get entangled, due to inter- 
actions. A measure of their entanglement at time t is 
given by the von Neumann entropy S(t) = — Tr(pi lnpi) 
where p\ is the reduced density matrix of particle 1 . The 
this way computed time dependence of entanglement is 
shown in Fig. [2]for three different values of h. A first ob- 
servation is that integrable and chaotic cases can hardly 
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tropics S c i (5, t) (with 5 = h) were computed, using a fixed 
partition of the 2-dimcnsional phase space in square cells 
of side y/S: 



FIG. 1: Surfaces of section for the Toda model at energy 
E — 7. (a) integrable case with masses mi = m,2 = 1; (b) 
mixed phase space for mi = 1, m,2 = 0.54. 



be told from each other by just comparing the degrees of 
entanglement to which they give rise. In this respect it 
is worth noting that the dependence of entanglement on 
initial conditions, which was observed in several papers 
{e.g. [1(3]) as the initial state was moved from the chaotic 
sea to an integrable island, is not a signature of chaos. 
Instead it is merely due to a decrease of the number of 
accessible quantum states, as one approaches the peri- 
odic orbit in the center of the island. Indeed, a similar, 
significant decrease of entanglement is recorded even in 
the completely integrable case of our model (not shown), 
starting with initial conditions that approach the central 
periodic orbit of the stable island; e.g., along a vertical 
line (qi — 1.2) in Fig. la. 

However, the most interesting observations are: first, 
as h decreases, the curves S(t, H), which describe the en- 
tanglement of the quantum Toda system as a function of 
time t at a given value of h, come closer and closer to a 
nontrivial limit curve (Fig. [2j full curve). Second, this 
curve is well approximated by a purely classical construc- 
tion, to be described next. 

From the classical Gaussian ensemble in the 4- 
dimensional phase space that corresponds to the initial 
quantum coherent state we have generated a sample of 
M initial points. To compute the evolution of this initial 
ensemble we have integrated the classical Hamilton equa- 
tion, using a fourth-order Runge-Kutta routine. The 4- 
dimensional ensembles thus computed at times t > were 
projected onto "reduced" ensembles in the (qi,pi) plane. 
To obtain a classical analog for the Von Neumann entropy 
of quantum reduced states, the Boltzmann-Shannon en- 
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where Wi(t) is the number of points in the reduced en- 
semble , which were found in the i'th cell at time t. We 
used enough cells to guarantee hosting of all points at all 
times. This classical entropy is expected to somehow ap- 
proximate the entropy of the quantum reduced ensemble, 
provided the latter ensemble mixes a sufficiently large 
number of states. 

We numerically observe that, as S — > 0, the curve 
S cl (8,t) approaches a limit curve Soo(£)- This seemingly 
crude classical construction works surprisingly well; in- 
deed, the limit classical curve is more and more closely 
approached by the curves of quantum entanglement as 
h — > 0. The approximation is so good, that even finer de- 
tails of the quantum curves, such as superimposed almost 
regular oscillations (fig. 2), can be explained in terms of 
classical phase space structures, as we shall presently ex- 
plain. 

Before coming to that, it is important to stress that 
our findings do not involve any contradiction. Mea- 
sures of quantum entanglement need not tend to in the 
h — > limit, as shown by simple examples: e.g., a system 
of two identical harmonic oscillators, in the pure state 
-j^Ho <g> H + -j^Hi <g> Hi, where H , Hi are the ground 
state and the 1st excited state of a harmonic oscillator, 
has entanglement entropy equal to ln(2), independently 
of the value of h > 0, even though in the limit h —> its 
Wigncr function tends to a Dirac 5 supported in the ori- 
gin of the 4-dimensional phase space. On the other hand, 
the meaning of the above described classical construction 
is an interesting question in purely classical terms, that 
may be related to a recently proposed classical separabil- 
ity entropy fl3j | . 

In Figs. 3,4 we show the numerically computed classi- 
cal and quantum entropies respectively. All the shown 
curves share some general qualitative features: notably, 
at long times they approach saturation values, via an 
average monotonic trend with strong oscillations super- 
imposed. Such oscillations look surprisingly regular even 
in chaotic cases, and gradually decrease in time; in the 
classical cases, at saturation they are on the order of 
fluctuations due to finite statistics. As the cell area 6 is 
decreased, for times less than a time scale to{5) a limit 
curve is approached, which is different in the regular and 
in the chaotic cases. The saturation time to may be 
thought of as the time when a finite 6 curve departs from 
the limit behaviour. For times t > to(S), entropy curves 
approach a saturation value, which scales proportional 
to l/\/S in the regular case and to ln(l/<5) in the chaotic 
case (see insets in Fig. 3). To infer estimates of saturation 
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FIG. 2: Entanglement entropy vs time for E — 7 and decreas- 
ing values of ft in a regular (a) and in a chaotic (b) case. The 
full curve is classical entropy. 



times from numerical data, tr> {$) was defined as the time 
when the average growth (represented by dashed lines) of 
the limit curve reaches the saturation value at the given 
5. 

Such numerical observations are easily explained by inter- 
preting entropy as the logarithm of an effective number 
of populated cells. To this end, one may first replace the 
initial Gaussian ensemble by a single 4-dimensional cell 
of size Vo. Saturation is then an obvious consequence of 
the fact that only a finite number of cells are available, 
due to energy conservation and compactness of the en- 
ergy shell. The saturation value of entropy is expected to 
scale with 6 like the logarithm of that number, i.e., in the 
leading order, like Cln(l/5) with (7=1, both in regular, 
and in chaotic cases. This is consistent with numerically 
observed values C ~ 0.86 and C = 0.93 for the regular 
and the chaotic case respectively. 

The average growth of entropy at large times, yet still 
far from saturation, may be understood similarly. The 
initial cell intersects different 3-dimensional manifolds of 
constant energy. In the chaotic case, evolution stretches 
the 3-dimensional section at energy E along the unsta- 
ble direction and contracts it along a stable one , at an 
exponential rate given by the maximal Lyapunov expo- 
nent Xe- Hence the number of cells which are popu- 
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FIG. 3: Classical entropy S c i(t) vs time for different values 
of the cell area S for E — 7. Above, an integrable case; the 
dashed line is given by const + 21n(t). Below, a chaotic case; 
the dashed line is given by const + t/10. Insets: saturation 
time to (circles) and saturation value S c i (squares) versus S. 
In the regular case, to ~ const + l/V~8, S c i ~ const + ln( 1/8). 
In the chaotic case: tz> « const + 101n(l/<5); S c i ~ const + 
ln(l/<5). 



lated by the projection of this section grows roughly like 
cxp(Aet). The number of cells populated by the projec- 
tions of all sections grows like a continuous superposition 
of exponentials with rates Xe varying in a small inter- 
val of size oc so, in the leading order, entropy grows 
proportional to t times the average Lyapunov exponent, 
independently of 8. Equating this to the saturation value 
the order of magnitude estimate to ~ A _1 ln(l/<5) is ob- 
tained, consistently with numerical results. 

In the regular case, the initial 4-dimensional cell inter- 
sects several invariant 2-tori. Each 2-dim section glides 
uniformly in time along the torus where it belongs, car- 
ried by the linear flow which is associated with the fre- 
quency vector of that torus; so its projection on the single 
particle phase space is just deformed quasi-periodically 
in time, with no average growth of the number of pop- 
ulated 2-dim cells. The number of 2-dimensional cells 
populated by all such projections nevertheless increases, 
because different sections move with different frequency 
vectors due to nonlinearity; so, although all sections were 
originally piled up over the same 2-dim cell, their projec- 
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FIG. 4: Entanglement entropy S(t) vs time for different val- 
ues of h, in an integrable case (above) and in a chaotic case 
(below). The continuous curves are the classical entropy S c i 
for S = 0.02 as in fig 4. Insets: saturation time to (circles) 
and saturation value S (squares) of entanglement, versus h. In 
the regular case: to ~ const + 1/VH, S « const + In (1/h). In 
the chaotic case: to ~ const + 10 ln(l/ft); S ss const + ln(l/ ft) 

tions slide apart linearly in time, oc ty/S (because the 
mismatch of frequency vectors is on the order of the lin- 
ear size of the cell) . As sections are a 2-parameters family 
, the 2-dimensional area covered by such projections in- 
creases quadratically in time : ^const. x St 2 in the leading 
order. Hence, the number of occupied cells of area S in- 
creases ~const.xi 2 , and entropy grows , in the leading 
order, like 21n(i), independently of 5. At time tn this 
should be on the order of magnitude of the saturation 
value ln(l/<5), whence the observed scaling of the satura- 
tion time to ~ 1/VS follows . 

As shown in Fig. 4, the behavior of quantum entangle- 
ment is quite similar. The same scaling with h of devia- 
tion times and saturation values emerges as the classical 
limiting curves are approached, both in the integrable 
and in the chaotic case fill ] - 

The oscillatory pattern which is seen in all figures, su- 
perimposed on the just discussed average behavior, is 
qualitatively explained as follows. In the regular case: 
the projection of each invariant torus on the (<?i,pi) plane 
has caustics at its boundaries, because there the tangent 
plane of the torus is "vertical" . As sections of the initial 



ensemble slide around their tori, the areas of their pro- 
jections change quasi-pcriodically in time, and in par- 
ticular become quite small anytime they come close to 
caustics. So long as different sections of the initial en- 
semble have not yet strongly dephascd from one another, 
this produces a regular sequence of minima in the num- 
ber of occupied cells. In the chaotic case: observation 
of the projected ensembles in the (qi,pi) plane, at times 
which correspond to minima in the regular sequence of 
oscillations, reveals that phase-space points are then con- 
centrated near a caustic-like line, that betrays presence 
of remnants of a broken torus. The projected ensem- 
ble appears to quasi-periodically concentrate near this 
line, until the 4-dimensional ensemble eventually leaks 
through the broken torus to invade most of the energy 
shell. 

We finally note that, while the quantum reduced entropy 
is the same for both particles, in general our classical 
reduced entropies approach each other only in the limit 
5->0. 

In conclusion, we have shown that there is a close similar- 
ity between the laws which rule the dynamical generation 
of quantum and classical reduced entropies. The classi- 
cal laws admit of simple heuristic explanations, which 
provide clues to understanding the quantum ones. Such 
arguments predict differences between the ways entan- 
glement is generated in classically regular and chaotic 
regimes, that should be observable for generic quantum 
multi-partite systems in the quasi-classical regime. 
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